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Abstract 

We present the Suzaku results of the mixed- morphology supernova remnant W28. The X-ray spectra 
of the central region of W28 exhibit many bright emission lines from highly ionized atoms. An optically 
thin thermal plasma in collisional ionization equilibrium, either of single-temperature or multi-temperature 
failed to reproduce the data with line-like and bump-like residuals at the Si Lymana energy and at 2.4- 
5.0 keV, respectively. The bumps probably correspond to radiative recombination continua from He-like 
Si and S. A simple recombining plasma model nicely fit the bump structures, but failed to fit low energy 
bands. The overall spectra can be fit with a multi-ionization temperature plasma with a common electron 
temperature. The multi-ionization temperatures are interpreted as elemental difference of ionization and 
recombination timescales. These results prefer the rarefaction scenario for the origin of the recombining 
plasma. 
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1. Introduction 

The evolution of thermal plasmas in supernova rem- 
nants (SNRs) is a dynamical time-sequence of shock- 
heated gas. At first, a large fraction of the expanding 
velocity is converted to random velocity (thermal energy). 
The shock-heated temperature is proportional to the mass 
of the plasma particles (electrons and ions); the ion tem- 
perature is more than 1000 times higher than the electron 
temperature (hereafter, Tg). Subsequently, the energy of 
ions is transferred to the electrons, and hence Tq gradu- 
ally increases. Then the high temperature electrons slowly 
ionize the ions to higher ionization states. Here we define 
the ionization temperature (Tz) as a parameter to repre- 
sent populations of ionization states; Tz is equal to the 
plasma temperature in collisional ionization equilibrium 
(CIE), which realizes the observed populations of ioniza- 
tion states. 

The typical timescale of ionization in SNR plasmas 
(density is ~ 1 cm~^) is ~ 10"^^"-^^ s (Smith & Hughes 
2010). Accordingly, most of the young-intermediate SNRs 
have lower Tz than Te, which is called an under- ionized or 
ionizing plasma, because ionization dominates over recom- 
bination. In fact, young-middle- aged shell-like SNRs have 
been generally described with an ionizing plasma model 
(e.g. Kinugasa & Tsunemi 1999). 

Recently, Suzaku has discovered strong radiative recom- 
bination continua (RRCs) in the X-ray spectra of three 
SNRs, IC443, W49B, and G 359.1-0.5 (Yamaguchi et al. 
2009; Ozawa et al. 2009; Ohnishi et al. 2011). The strong 
RRCs appear when recombination dominates over ioniza- 
tion (recombining plasma). In fact, the X-ray spectra of 
these SNRs can be fit by neither CIE plasma, ionizing 
plasma, nor any combination of them. Detailed spectro- 



scopic analysis with the recombining plasma code revealed 
that Tz are 1.5-3 times higher than Tg. These SNRs have 
some common features. (1) All are categorized into mixed- 
morphology SNRs (MM SNRs: Rho & Petre 1998), which 
have a radio shell with centrally- filled thermal X-rays. 
(2) These are associated with shocked molecular clouds 
and GeV/TeV 7-ray emissions (e.g. Yusef-Zadeh et al. 
2003a; Abdo et al. 2010a). 

W28 (G 6.4-0.1), a bright MM SNR, exhibits above 
two conditions; a shock interaction with the ambient gas 
indicated by the expanding Hi shell (Velazquez et al. 
2002), broad CO, CS, and H2 emission lines (Arikawa et 
al. 1999; Reach et al. 2005), and OH (1720 MHz) masers 
(Yusef-Zadeh et al. 2003b). Also GeV and TeV 7-rays 
were detected (Abdo et al. 2010b; Aharonian et al. 2008). 

From the estimated distance of 1.9±0.3 kpc (Velazquez 
et al. 2002), the angular size of ^50^ is ~30 pc. The 
large size, possible association to the nearby pulsar (Kaspi 
et al. 1993), and dynamical evolution of the Hi shell 
(Velazquez et al. 2002) constrain the age of W 28 to be 
33,000-150,000 yr. 

ROSAT and ASCA mapped the entire remnant in the 
X-ray band (Rho & Borkowski 2002) and found center- 
filled structure with partial shells at northeast and south- 
west. The X-ray spectrum in the central region had hard 
X-ray emissions up to ~ 7 keV with the iron (Fe) Ka line, 
and was fit with a two-temperature plasma of 0.67 keV 
and 1.8 keV in CIE or ionizing. Thus no hint of recom- 
bining plasma has been reported from W 28, although this 
SNR shares the common features with the three SNRs of 
recombining plasmas. The aim of this paper is to search 
for and quantitatively study recombining plasmas with 
the superior energy resolution and the large effective area 
of the X-ray Imaging Spectrometer (XIS: Koyama et al. 
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Table 1. Log of Suzaku observations of W28. 



Sequence no. 


Aim point 


Start date 


Effective 


Field name 




a (J2000.0) S (J2000.0) 




exposure 




505005010 


18^00"^17" -23°21'59'' 


2010/04/03 


73.0 ks 


Center 


500008010 


18^03"^49" -22°01W 


2006/04/07 


40.7 ks 


Nearby sky 




Galactic longitude (degree) 



Fig. 1. Wide- field X-ray image retrieved from the ROSAT All Sky Survey in grayscale (Snowden et al. 1995). Overlaid contours 
are the radio continuum map at 1.4 GHz by the Very Large Array (Dubner et al. 2000). The XIS FOVs are shown with the solid 
squares. 



2007) onboard Suzaku (Mitsuda et al. 2007). Throughout 
this paper, the distance of 2 kpc is adopted and statistical 
errors are at the 90% confidence level. 

2. Observations and Data Reductions 

We observed the central bright region of W 28 with the 
XIS as a part of the recombining plasma survey project 
(PI: Koyama, K.). The local background was obtained 
from the Suzaku archive near W28. The locations of the 
two fields and the observation log are respectively shown 
in figure 1 and in table 1. 

The XIS consists of four CCD cameras each placed 
at the focal planes of four X-Ray Telescopes (XRTs: 
Serlemitsos et al. 2007). Three sensors employ Front- 
Illuminated (FI) CCDs (XISO, 2, and 3), while the other 
employs a Back-Illuminated (BI) CCD (XIS 1). The entire 
region of XIS 2 and one of the edge region of XIS have 
not been functional since anomalies in 2006 November and 
in 2009 June, respectively. The field of view (FOV) of the 
XIS combined with XRT covers a 18' x 18' region with 
the pixel scale of 1" pixel" ^. The angular resolution of 
l.'9-2.'3 in the half-power diameter is almost independent 
of photon energies and off- axis angles within ~ 10%. The 
total effective area of the operational XIS 0, 1, and 3 com- 
bined with the three XRTs is 1070 cm^ at 1.5 keV. Due 
to the low orbital altitude of Suzaku at ~550 km, the XIS 
achieves a low and stable background environment. 

The XIS observation was made with the normal clocking 
mode. The software package HEASoft version 6.11 and 



the pipeline processing version 2.4 were used for the data 
reduction. To restore the radiation-induced degradation 
in the energy gain and resolution, the spaced-row charge 
injection technique (SCI: Bautz et al. 2004) was applied 
with the makepi files version 20110621 provided by the 
XIS team (Uchiyama et al. 2009). Then the systematic 
uncertainty in the energy scale is ^ 10 eV at 5.9 keV. 
We removed hot and flickering pixels, and events during 
the South Atlantic Anomaly passages and in the Earth 
night-time and day-time elevation angles below 5° and 
20°, respectively. 

3. Analyses and Results 

3.1. Background Estimation 

Figure 2 shows the X-ray images of the W 28 center and 
the nearby-sky fields. The non-X-ray background (NXB) 
generated by xisnxbgen (Tawa et al. 2008) was subtracted, 
and then the vignetting and exposure corrections were 
made. 

We can see excess emission from W28 in the soft X-ray 
band, but no excess in the hard band (figure 2). We ex- 
tracted the source and background spectra from the solid 
and dashed circles in figure 2, respectively. The results 
are shown in figure 3, where the effective area of the back- 
ground is normalized to that of the source region. 

Compared with the background, the source spectrum 
shows obvious excess below ^ 5 keV but no excess in the 
hard X-ray band ^ 5 keV (figure 3), consistent with the 
images. This result is, however, inconsistent with the pre- 
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(a) Soft 




8.2 8.0 6.6 6.4 



Galactic longitude (degree) 

(b) Hard 




8.2 8.0 6.6 6.4 



Galactic longitude (degree) 

Fig. 2. Band-limited images with the XIS: (a) soft X-rays in 
0.5-5.0 keV and (b) hard X-rays in 5.0-8.0 keV. The data 
with the three CCDs were merged. The NXB was subtracted 
and then the exposure and vignetting effects were corrected. 
The source and background extraction regions are shown in 
the sohd and dashed circles, respectively. 



Si 




Energy (keV) 

Fig. 3. Comparison of the source (blue) and the normalized 
background (red) spectra. For visibility, only merged spec- 
tra of FI (XIS and XIS 3) are displayed. The solid and 
dashed black lines show the center energies of the Ka lines 
from He-like and H-like ions, respectively. The gray line indi- 
cates the CXB spectrum. 

vious ASCA observations, which claimed the detection of 
hard X-ray emission up to ~ 7 keV with the Fe Ka line 
(Rho & Borkowski 2002). 

One problem is position-to-position fluctuations of the 
cosmic X-ray background (CXB). We therefore estimated 



the CXB flux assuming the photon index, the average sur- 
face brightness, and the Galactic absorption to be 1.412, 
6.38x10"^ erg cm-^ s"^ sr"^ (Kushino et al. 2002), and 
1.4 X 10^^ H cm~^ (Dickey & Lockman 1990), respectively. 
Then the expected CXB level (the gray curve in figure 3) 
is only about 5% of the total background flux. Thus the 
error due to the CXB fluctuation can be neglected. 

The major source of the background is the Galactic dif- 
fuse X-ray emission (GDXE) because W 28 is located near 
the Galactic plane and the Galactic center (GC). In fact, 
the Fe Ka line at 6.7 keV seen in figure 3 is due to the 
GDXE. Since the surface brightness of the GDXE quickly 
decreases with the increasing distance from the Galactic 
plane with the scale height of ^ 1° (Uchiyama 2010^), 
we selected the background region from the same Galactic 
latitude (figure 2 left) of W28 (figure 2 right). The GDXE 
uncertainty is, therefore, due to the longitude distribution 
along the plane. The flux of the GDXE decreases as the 
distance from the GC increases with the scale length of 
-30° (Uchiyama 2010). Then the GDXE on the back- 
ground region is estimated to be about 5% smaller than 
that on W28. Still we found no excess from W28 above 
the background in the hard X-ray band, and hence no hard 
X-ray emission above 5 keV for W 28 is very conservative 
conclusion. 

3. 2. Spectrum 

We subtracted the background spectrum (the red data 
in figure 3) from the source spectrum (the blue data in fig- 
ure 3). The result is shown in figure 4. We can see many 
K-shell emission lines from highly ionized neon (Ne), mag- 
nesium (Mg), silicon (Si), and sulfur (S) as well as L-shell 
emission lines from Fe. We therefore examined several 
optically thin thermal plasma models. In the spectral 
analysis, we used the SPEX (Kaastra et al. 1996) soft- 
ware version 2.02.02. The redistribution matrix and the 
auxiliary response functions were generated by xisrmfgen 
and xissimarfgen (Ishisaki et al. 2007), respectively, with 
the latest rmfparam files version 20110708. The abun- 
dances for the absorption of interstellar medium (ISM: 
Morrison & McCammon 1983) were assumed to be one 
solar (Anders & Grevesse 1989), using the Absm model 
in SPEX. The abundances for Ne, Mg, Si, S, and Fe in 
the thermal plasma were free parameters, while those of 
oxygen (O), argon (Ar) and calcium (Ca), and nickel (Ni) 
were assumed to be the same as those of Ne, S, and Fe, 
respectively. The other elements, including carbon (C) 
and nitrogen (N) were fixed to one solar. 
3.2.1. CIE Plasmas 

We tried a single-temperature optically thin thermal 
plasma model in CIE (the Cie model in SPEX). The re- 
sult is shown in figure 4a. Although we fit simultaneously 
for the FI and BI spectra, we show only the best-fit FI 
spectrum for simplicity. As is seen in figure 4a, this model 
obviously failed to reproduce the spectrum with x^/d.o.f. 
of 2160/641. We therefore examined a two- temperature 
CIE model (figure 4b), with coupled abundances for each 

1 Available at http://repository.tksc.jaxa.jp/pl/dr/IS8000028000/en 
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Fig. 4. Background-subtracted spectra (the black crosses) with the best-fit models. For visibility, only data and models for the 
merged FI are displayed except for (f). Each panel shows results of different models: (a) single-temperature CIE, (b) two-temperature 
CIE, (c) single- Tz recombining plasma model, (d) multi-Tz recombining plasma model, (e) the same as (d) but artificial Gaussians 
at 0.8 and 1.2 keV (pink) are added, and (f) the same as (e) but for BI. The thick gray curves are the total emission, while the thin 
colored curves are each component: low-A^T component in dark red and high-/cT component in dark green for (b); C and N in red, 
O and Ne in orange, Mg in light green. Si in light blue, S, Ar, and Ca in dark blue, Fe and Ni in purple for (d), (e), and (f). The 
solid arrows in (c) indicate the Ka line centers of He-like Mg and Ne, and L-shell complex of Fe. The solid arrows in (d) are L-shell 
transitions of Fe. 



element between the two CIE components. This model 
also failed to reproduce the spectrum with x^/d.o.f. of 
1753/639. The best-fit parameters for these CIE mod- 
els are given in table 2. In addition to the 0.77-keV 
plasma, W28 requires 0.24-keV plasma instead of the 1.8- 
keV plasma of the ASCA results (Rho & Borkowski 2002). 
The authors proposed that W 28 is a pecuhar MM SNR 
because of the presence of the high-temperature plasma 
of 1.8 keV. Our results of the 0.24+0.77-keV tempera- 
ture plasma revised the interpretation of W28 to be a 
typical MM SNR (0.4-0. 9-keV temperature; Rho & Petre 



1998). This revised picture comes from the absence of the 
high energy photons above ~ 5 keV in the Suzaku spec- 
trum, which would be due to proper background subtrac- 
tion (subsection 3.1). To improve the value of x^/d.o.f., 
we added further CIE components, but found no signifi- 
cant improvement. The most remarkable residuals in any 
CIE model are line-like excess at the Si Lymana energy 
(2.0 keV) and bump-like structures at 2.4-5.0 keV. 
3.2.2. Recombining Plasma Model with Single 

The bumps at 2.4-5.0 keV have two rising edges at 
2.4 and 3.2 keV, which are most probably due to RRCs 
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Table 2. Best-fit parameters with CIE and single- Tz recombining plasma models. 



Parameter 



Single-temperature CIE Two-temperature CIE Single-Tz recombining plasma 



A^H (lO^i cm-2) 
VEMi (10^^ cm- 
VEM2 (10^^ cm" 
kTi (keV) 
kT2 (keV) 
kTe (keV) 
kT^ (keV) 



Q Q0+U.U7 

o.yu_Q Q7 

21 1+^-^ 
^-'-•-'--0.6 



62+°-°^ 
'-'•'^^-0.01 



6.28"!^ 
103l{^ 

13.61^0:^5 
0.24+^-01 



+0.U4 
-0.03 
+0.9 
-0.9 



-3y 
-3y 



0.77 



0.01 
+0.01 
-0.01 



2.99 
64.6 



0.43 (fixed) 
0.95 (fixed) 



Abundance (solar) 
Ne 
Mg 
Si 

s 

Fe 



QQ+O-02 

u.oy_o_o2 
S6+^-°2 

U.OO_Q Q2 

'-'•^^-0.01 

^+0.04 
-0.04 

'-'•-'-^-0.01 



0.50] 



n+o-oi 

U.iO_Q Q;L 

U.4i_0.02 

n 09+0.02 
'-'•'^^-0.02 
4S+0-04 

23+°-^^ 

U.ZO_Q Q-|^ 



<0.01 
0.18 



+0.01 
-0.01 

0.20iroJ 
0.26iro^2 



2160/641 



1753/639 



6914/642 



: J UeUpdV , where n©, Up, and V, are the electron and proton densities, 



* Volume emission measure at the distance of 2 kpc: VEM - 
and the emitting volume, respectively. 



of He- like Si and S, respectively (Ohnishi et al. 2011). 
Together with the strong Si Lymana residual, the spec- 
trum is due to a recombining plasma, at least for Si and 
S. We thus tried a recombining plasma model in the RRC 
band of Si and S (1.7-5.0 keV). Then the residuals found 
in the 1.7-5.0 keV band disappear with the best-fit of 
0-95to:o9 keV and Te of OAStom keV. Fixing these best- 
fit temperatures, we made a full energy band fit. The free 
parameters were abundances of Ne, Mg, Si, S, and Fe, and 
A/'h- The result is, however, unacceptable with x^/d.o.f. 
of 6914/642 (figure 4c). As shown by the solid arrows 
in figure 4c, remarkable excess are found at the Ka lines 
of He-like Mg (1.35 keV) and Ne (0.92 keV) and the Fe 
L-shell complex at '^0.8 keV. 
3.2.3. Multi-Tz Recombining Plasma 

The data excess near the Ka lines of He-like Mg and Ne 
indicates underestimation of the He-like line fluxes com- 
pared to the H-like line (see the solid arrows in figure 4c). 
This supports that Tz for Mg is significantly lower than 
those of Si and S. We therefore introduced a recombin- 
ing plasma model with different Tz for each element (fig- 
ure 4d). Since no spectral code to directly describe such 
a plasma is available, we approximated the model with 
the combination of several single-Tz plasmas whose Tq are 
coupled among the all plasmas. We grouped the elements 
into six sub-groups assuming the same abundances in each 
sub-group (see the caption of figure 4). For C and N, we 
fixed Tz to Tq. 

As shown in figure 4d, this model (model A) well repro- 
duced the overall structure of the observed spectrum. Still 
some residuals remain at ~ 0.8 and ~ 1.2 keV (see the solid 
arrows in figure 4d). These features are known to origi- 
nate from incomplete atomic data in the current spectral 
code. The former comes from the uncertain intensity ratio 
of the Fe L-shell transitions of 3s^2p and 3d^2p (Gu et 
al. 2007 and reference therein). The present code may un- 
derestimate the 0.8 keV line flux. The latter is due to the 



Table 3. Best-fit parameters with a multi-Tz recombining 
plasma model. 



Parameter 




B 


A^H (10^^ cm-2) 
VEM (10^^ cm-3)* 
kTe (keV) 


4.65l-.i^2 
60.7+?-? 

0.39l^o;'o'2 


4 70+0-1^ 
4T.lt^-i 


kTz (keV) 
... Ne 
... Mg 
... Si 
... S 
... Fe 


0.66tH^ 
1 ni+o-07 

-'-•'-'-'--0.06 

^8+0-04 

'-'•^^-0.04 


46+0-05 

'-'•'^'^-0.06 

1 01 +0-08 
-'-•'-'-'--0.08 

0.96l°;l^ 


Abundance (solar) 

... Ne 

... Mg 

... Si 

... S 

... Fe 


16+0-01 

U.iO_Q Q2 

n 09+0.02 

'-'•'^^-0.02 

26+0-02 

U.ZD_Q Q3 

n QQ+0.08 

'-'•'^'^-0.08 

10+0-01 


n 9Q+0.03 
'-'•^'^-0.02 

44+0-04 
S4+0-03 
40+0-09 
11+0-01 

'-'•-'--'--0.01 


xVd.o.f. 


1268/636 


893/634 



* The value at the distance of 2 kpc. 



lack of a number of the Fe L-shell lines by the transitions 
from highly excited states of n > 5 (Brickhouse et al. 2000). 
Thus, we artificially added two Gaussians at 0.8 keV and 
1.2 keV (figure 4e and 4f). Then the fit largely improved 
(model B; xV^-oi. = 893/634). The best-fit parame- 
ters for Model A and B are listed in table 3. We see that 
the best-fit Tz are essentially the same between model A 
and B. The best-fit Tz are different among elements with 
higher temperature than Tq (table 3). 
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Fig. 5. Time evolution of Tz in recombining plasmas with the initial conditions of (a) Tzo ■ 

realized by an ionizing plasma with = 10 keV and rietioni = lO"*^"*^ s cm~^. The gray circles with bars show the observed Tz with 
1-cr statistical errors. The direction of evolution as a function of Uetrec is indicated by the vertical arrows. The solid lines in red and 
blue respectively show the initial and final states, while the dashed blue lines show transitional states. 



4. Discussion 

We have found that the observed spectrum can be re- 
produced by a multi-Tz recombining plasma. The recom- 
bining plasmas in SNRs are not expected in the standard 
scenario, and hence there must be a missing branch of evo- 
lution of SNRs. Several ideas have been proposed, such 
as rapid electron cooling by adiabatic rarefaction (Itoh & 
Masai 1989), or conduction cooling by ambient cold mat- 
ter (Kawasaki et al. 2002), extra- ionization by high-energy 
photons (Kawasaki et al. 2002) or by supra-thermal par- 
ticles (Ohnishi et al. 2011), and so on. Any scenario may 
compose of two processes; (1) initial process to make a 
plasma of higher Tz than Tq , and (2) gradual relaxation 
of the recombining plasma toward CIE. We propose that 
the process (2) carries the essential roll to make element- 
dependent Tz. Based on the simulation of the Tz relax- 
ation in time and comparison with the observed multi- 
Tz values, possible origin of the recombining plasma from 
W28 is discussed in the following subsections. 

4.1. Recombination Timescale 

Difference of Tz among elements can appear in a relax- 
ation epoch of recombining plasma to CIE because the 
relaxation timescale is different from element to element. 
Numerical studies have demonstrated that the relaxation 
timescale basically shows monotonic increase with atomic 
number, except for Fe and Ni (e.g., see figure 1 of Smith 
& Hughes 2010). These two heavy elements have shorter 
timescales than the other lighter elements at Tz from sub 
to a few keV, because the L-shell processes are dominant 
for Fe and Ni, while the K-shell processes are dominant 
for the others in this temperature range. The observed Tz 
values appear to satisfy this trend (table 3). 

To examine whether the difference in recombination 
timescale can quantitatively explain the observed Tz, we 
simulate time evolutions of Tz in a recombining plasma 



by using the Neij model in SPEX. The parameters are the 
initial ionization temperature Tzo, Tg, and the recombi- 
nation time Uetrec^ whcrc He and tree are the density of 
electrons and the elapsed time, respectively. 

We have tried many sets of (Tzo, tree) ^oi fixed Te of 
0.4 keV following the best-fit value (table 3). We obtained 
similar best-fit Tz distribution from different parameter 
set of (Tzo, tree)- We therefore constrain the parame- 
ter set within physically reasonable values. One plausible 
case of Tzo = 3 keV is given in figure 5a, which shows the 
simulated decrease of Tz with time (netree)- We find the 
best-fit recombination time to be Uetrec = 10^^*^ s cm~^. 
The simulated Tz values after 10^^'^ s cm~^, is close to 
the observation with x^/d.o.f. of ~ 21/4. We caution 
that the x^/d.o.f. value is just for reference to compare 
a more complex simulation given in the next subsection, 
and should not be taken seriously, because the major un- 
certainty of the observed Tz is not statistic error but sys- 
tematic one due to incomplete modehng and/or atomic 
data in the analyses. 

4.2. Initial Variation of Tz 

In a more general case, initial Tzo should be different 
among elements. We hence simulate the time evolution 
of Tz in this case. The simulation consists of two steps. 
In the first step, we simulate the evolution of Tz in an 
ionizing plasma. This gives the element-dependent Tzo- 
Then using this result, we simulate a recombining plasma 
as was done in the previous subsection. 

The parameters are the electron temperature during the 
ionizing phase T^, the ionization time netioni, and the re- 
combination time netree- The electron temperatures dur- 
ing the two processes, T^ and Tg, are assumed to be 10 keV 
and 0.4 keV, respectively. 

In figure 5b, we plot decrease of Tz during the recom- 
bining phase for the case of T^ = 10 keV and netioni = 
10^^"^ s cm~^. This combined process gives a better fit 
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to the observed T^, than that of the single relaxation 
process of recombining plasma (see subsection 4.1) at 
rietrec = lO^^'^ s cm"^ with xVd.o.f. of - 17/4. One 
may argue that other combinations of the parameter val- 
ues may also work as well. As we already noted, many 
combinations were examined and the recombination times 
are always found to be ^10^^-^ s cm~^ or longer. Thus, 
the result of Uetrec does not depend much on the initial 
conditions of T^o- 

4.3. Origin of Recombining Plasma in W 28 

We have shown that the combination of the processes 
(1) and (2) provides a plausible scenario for producing the 
observed spectrum with the initial conditions given in fig- 
ure 5(b) (see the beginning of Discussion and section 4.1). 
The initial conditions can be set by collisional ionization 
in thermal pool of electrons with no special ionization pro- 
cess such as photo-ionization or supra-thermal electrons. 
This supports the electron cooling scenario for W28. In 
the simulation, ne^rec is ~ 10^^'^, which corresponds to 

10^ X (ne/1 cm~^)~^ yr, roughly consistent with the age 
of W28. Thus, the recombination time is consistent with 
the rarefaction scenario (Itoh & Masai 1989). 

This recombination time is also nearly the same order 
of the cooling timescale by thermal conduction, which is 
estimated to be ~2xl0^ x (ne/lcm~^) yr for the plasma 
of 6-pc diameter (see equation 5 in Kawasaki et al. 2002). 
However, we found no spatial correlation of the recombin- 
ing plasma to the molecular clouds (Sawada et al. 2012, 
in preparation) , which may not favor the thermal conduc- 
tion scenario. We therefore propose the rarefaction sce- 
nario (Itoh & Masai 1989), and discuss along this scenario 
in the following paragraphs. More quantitative compari- 
son of these two models will be given in the next paper 
using the results of spatial distribution of the recombining 
plasma (Sawada et al. 2012, in preparation). 

The initial event of the rarefaction scenario would 
be a type-II supernova explosion in dense circum-stellar 
medium (CSM) made by stellar winds from the massive 
progenitor. Assuming a spherical plasma with the diame- 
ter of 6 pc, the mass of the X-ray emitting gas in the center 
region is estimated to be ~5 Mq. This value is similar to 
that of a massive progenitor, and hence provides further 
circumstantial support for the rarefaction scenario. The 
shock-heated ejecta and CSM experience rapid electron 
cooling due to adiabatic expansion after shock break-out 
to the rarefied ambient ISM. The break-out and hence 
the cooling, occur at the early phase of the SNR evolu- 
tion (about several 10-100 yr after the explosion; Itoh 
& Masai 1989). Thus most of the life would be in the re- 
laxation phase of the recombining plasma, consistent with 
the large rig tree values obtained in our simulations. Due to 
the dense environment of CSM, the shocked gas quickly 
reaches a temperature of several keV before the break- 
out, in good agreement with the Tzo values assumed in 
our simulations. 

Shimizu, Masai, & Koyama (2011) extended the previ- 
ous calculation of Itoh & Masai (1989) to non-symmetric 
cases and found that center-filled X-ray morphologies are 



realized with anisotropic CSM due to the rotation of a 
progenitor star. The size of the X-ray emitting plasma is 
about a few 10 pc, which is consistent with that of W28. 
Therefore, this scenario can explain both the spectral and 
spatial structures of the X-ray emission from W28. 
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